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Summation of divergent power series by means of factorial series 
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Abstract 

Factorial series played a major role in Stirling's classic book Methodus Dijferentialis (1730), but now only a few 
specialists still use them. This article wants to show that this neglect is unjustified, and that factorial series are useful 
numerical tools for the summation of divergent (inverse) power series. This is documented by summing the divergent 
asymptotic expansion for the exponential integral Ei [z) and the factorially divergent Rayleigh-Schrodinger perturba- 
tion expansion for the quartic anharmonic oscillator. Stirling numbers play a key role since they occur as coefficients 
in expansions of an inverse power in terms of inverse Pochhammer symbols and vice versa. It is shown that the 
relationships involving Stirling numbers are special cases of more general orthogonal and triangular transformations. 

Keywords: Factorial series; Divergent asymptotic (inverse) power series; Stieltjes series; Quartic anharmonic 
oscillator; Stirling numbers; General orthogonal and triangular transformations; 
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1. Introduction 



Power series are extremely important analytical tools not only in mathematics, but also in the mathematical treat- 
ment of scientific and engineering problems. Unfortunately, a power series representation for a given function is 
\jQ ' from a numerical point of view a mixed blessing. A power series converges within its circle of convergence and 

\^ , diverges outside. Circles of convergence normally have finite radii, but there are many series expansions of consider- 

■^f ' able practical relevance, for example asymptotic expansions for special functions or quantum mechanical perturbation 

expansions, whose circles of convergence shrink to a single point. 
lO ' The summation of divergent (inverse) power series is an old problem of mathematics, which is of considerable rel- 

^P , evance also in related disciplines. Many different summation techniques have been developed which are often capable 

of associating a finite value to a divergent series. A highly condensed overview of various summation techniques was 
recently given in [i54i. A ppendices A and B]. The role of summation techniques in theoretical physics was discussed 
in the recent review 11711 . 

The topic of this article is the summation of divergent power series via so-called /acfona/ series. A factorial series 
;h ' for a function Q.: C ^- C, which vanishes as z — > 4-°°, is an expansion of the following type: 



^" z z(z + l) z(z+l)(z + 2)^'" ioWv+r '' 

Here, (z)v+i = r(z + V + l)/r(z) = z{z + 1) . . . (z + v) is a Pochhammer symbol. The function £l{z) represented 
by the factorial series dl.lb may have simple poles at z = —m with m G Nq. The definition (II. lb is typical of the 
mathematical literature. It will become clear later that the separation of the series coefficients into a factorial «! and a 
reduced coefficient a„ often offers formal advantages. 

The use of factorial series for the summation of divergent (inverse) power series is not a new idea. There is 
Watson's classic article on the transformation of an asymptotic inverse power series to a convergent factorial series 
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1 4511 . This topic was also considered in articles by Nevanlinna BOll . Malgrange ||28|] . and Ramis and Thomann |13( 
Thomann |42, 43] discussed the summation of formal power series with the help of generalizations of factorial series, 
where the Pochhammer symbols (z)„+i — z{z+ I) . . . (z + n) are replaced by products z(z + fi) ...{z + t„). There is 
also a recent article by Delabaere and Rasoamanana 02011 on the connection of Borel summation and factorial series, 
which can be viewed to be an extension of a short discussion in Borel's classic book |6, pp. 234 - 245]. 

These references show that there is no lack of knowledge about factorial series in general and about their use as 
summation tools. Unfortunately, this knowledge is restricted to a relatively small group of specialists, and there is a 
deplorable lack of public awareness about factorial series and their numerical usefulness. My claim is supported by the 
fact that factorial series are not mentioned in the 2010 Mathematics Subject Classification (MSC), in the web database 
MathWorld (http : //mathworld . wolfram . com), or in Wikipedia (http : //en . wikipedia . org/wiki/Main_Page). 
This neglect is not justified, and I believe that the full potential of factorial series as summation tools has not yet been 
reaUzed. 

Section 12] describes how I had become interested in factorial series. In Section[3] the basic properties of factorial 
series are reviewed. In Section H] it is shown that inverse power series and factorial series can be transformed into 
each other with the help of Stirling numbers. In Section|5] the transformation of Stieltjes series is considered. Section 
|6]describes the transformation of the divergent asymptotic series for the exponential integral to a convergent factorial 
series. In Section|7l a transformation for power series, which is essentially a factorial series in 1 /z, is discussed, and in 
Section[8j this transformation formula is used for the summation of the divergent Rayleigh-Schrodinger perturbation 
expansions for the ground state energy eigenvalue of the quartic anharmonic oscillator This article is concluded by 
a short outlook in Section |9l Those properties of Stirling numbers, which are for our purposes most relevant, are 
reviewed in [Appendix A[ In [Appendix B| it is shown that the transformation formulas for inverse power and factorial 
series considered in Section[5] are just special cases of more general transformation formulas involving triangular and 
orthogonal matrices. 

2. My "rediscovery" of factorial series 

In 1985/1986 I became interested in Levin's sequence transformation |25], and I tried to understand the mathe- 
matical theory behind it. My interest was aroused by two articles by Smith and Ford |39, 40] who had shown that 
certain variants of Levin's sequence transformation were among the most powerful as well as most versatile sequence 
transformations known at that time. 



In my work on Levin's sequence transformation ||25[1 I discovered that its derivation becomes almost trivially 



simple if we start from the model sequence 0461 Eq. (3.2-9)] 

S„ = S + COnZn , n G No . (2.1) 

The remainder estimates a)„ are assumed to be known, and the correction terms z„ should be chosen in such a way 
that the products (X)„Zn provide sufficiently accurate and rapidly convergent approximations to the actual remainders 
{'"n}r=o of the sequence {in}"=o which is to be transformed. 

In this approach, only the correction terms {z„}'^^q have to be determined. If good remainder estimates can be 
found, the determination of z„ and the subsequent elimination of a)„z„ from s„ often leads to better results than the 
construction and subsequent elimination of other approximations to r„. 

The model sequence ( 12. Il l has another indisputable advantage: A sequence transformation, which is exact for this 
model sequence, can be constructed easily under very mild conditions. Let us assume that a linear operator T can 
be found which annihilates the correction term z„ for all n £ No according to T{zn) = 0. Then we obtain a sequence 
transformation, which is exact for the model sequence (12. Il l, by applying T to the ratio [s„ — s\/(i)n = Zn- Since T 
annihilates Zn and is by assumption Unear, the following sequence transformation £/' is exact for the model sequence 
(EBSEq. (3.2-11)]: 

^is„,co,,) = ^^=s. (2.2) 

r(i/(B„) 



The annihilation operator approach was introduced in 0461 Section 3.2] in connection with my rederivation of 



Levin's transformation [251, but I also constructed in this way some other, closely related sequence transformations 



kd . Sections 7-9]. Brezinski and Redivo Zaglia il5ill6ll and Brezinski and Matos il4|] showed later that this approach 



is actually much more general: The majority of the currently known sequence transformations can be derived via ( 12 .21 1 
(for further references on this topic, see [52, p. 1214]). 



We obtain a model sequence for Levin's sequence transformation 11251] by assuming that Zn in ( 12. Il l is a truncated 
power series in l/(j3 +«) iH, Eq. (7.1-1)]: 

^"^Ltr^' i3>0. (2.3) 

The A;th power of the finite difference operator A defined by Af{n) = f{n + 1) — /(n) annihilates an arbitrary poly- 
nomials Pk-i{n) of degree k—1 inn. Thus, the weighted difference operator t — A^{I5 +«)*^^ is the appropriate 
annihilation operator for z„ defined by ( 12.31 ). and Levin's sequence transformation 1,25,1 can in the notation of 1,46- Eqs. 
(7.1-6) and (7.1-7)] be expressed as follows: 

The numerator and denominator sums of this and of related transformations can also be computed recursively 115 2[ 
Section III]. 

This undeniable success inspired me to look for other applications of the annihilation operator approach. If we 
replace in (I2.3l l the powers (/3 + n)^ by Pochhammer symbols (j3 + n)j, we obtain a truncated factorial series in j3 + n: 

k-i c' 
^" = E7r^' i3>0. (2.5) 



Now, f = A*^(j3 +n)i-i is the appropriate annihilation operator, and we obtain fl46l Eqs. (8.2-6) and (8.2-7)]: 

y^ (p,s„,co„j — ,,,,„ ^^ ; — 7- ~ — , K,nGNo. (2.6) 

' An(^+«),-iM,] k ./A (^+„+j),_, 1 



My derivation of this sequence transformations, whose theory was developed in 11461 Section 8], was entirely 
based on heuristics. I had only looked for situations in which I could apply the annihilation operator formalism ( I2.2l l 
effectively. If the correction term z„ is according to (12.5b a truncated factorial series, it can be annihilated easily, but I 
had no idea whether the resulting sequence transformation ^^ ( j3 , i,, , ©„ ) would be computationally useful or not. 

o5^/_^(j3,s„, 0),,) was first used for the evaluation of auxiliary functions in molecular electronic structure calcula- 



tions H56I1 . Later, it was used with considerable success in the case of slowly convergent or divergent alternating series 
(numerous references are listed in 0521 p. 1225]). Currently, yi^"'{f},s„,CO„) is used quite a lot in optics i7l 4l3il26ll27ll . 



When I constructed y^" {P,Sn,co„) in 1986, I had no idea what the correction term ( 12.51 ) actually is: Factorial 
series had not been part of my mathematical training. It took a while until I found out that ( I2.5l l is a truncated factorial 
series, and that Nielsen's classic book [31], which is still one of the principal references on factorial series, had already 
been waiting for quite a while on my bookshelf. 

In my later work on convergence acceleration and summation processes, I noticed that ^1^" {P,s„,co„) and 

(n) 

yj^ (j3,i„,a)„) usually have similar, but not identical properties. Nevertheless, in some cases spectacular differ- 
ences were observed. For example, in summation calculations for the divergent Rayleigh-Schrodinger perturbation 



expansions of the ground state energies of anharmonic oscillators II55I1 . we observed that Levin's transformation 



.^" (/3 , s„ , co„ ) ultimately produces divergent results, whereas ^j- {p,s„,co„) produces very good results (see also 



kSIl ). A similar divergence of Levin's transformation was also observed by Cizek, Zamastil, and Skala il9tl . Needless 
to say that these observations puzzled me. 

The construction of Levin's transformation is based on the assumption that the ratio [s„ — s]/co„ can be expressed 
as an inverse power series, whereas ,S^^" {p,s„,co„) implicitly assumes that [s„ — s] /co„ can be expressed as a factorial 

series. For me, it was a plausible hypothesis that the observed differences of J^f^" {l3,s„,co„) and ,y^' {j5,s„,co„) 
could be related to different properties of inverse power and factorial series, respectively. Because of my complete 
lack of knowledge about factorial series, I first had to study their properties. The results presented in this article are in 
some sense a by-product of these studies. 

3. Basic properties of factorial series 

Factorial series have been known for a very long time. In Tweedle's annotated translations of Stirling's classic 
Methodus Differentialis it is remarked that Stirling was not the inventor of factorial series. Apparently, Stirling became 
aware of factorial series by the work of the French mathematician Nicole ll44l p. 174]. However, Stirling used factorial 
series extensively and thus did a lot to popularize them. 

The application of higher powers of the finite difference operator A = A^ to a factorial series in z yields an ex- 
tremely compact result. If we use A^[«!/(z)„+i] = { — lY {n + k)\ / {z)„+ii+i with k G No, we obtain 

^'Mz) = i4^^ (-1)^ E "-f^ = (-1)^ E ^. (3.1) 

v=0 \^)v+l v=0 {^)v+k+l f^^f, [Z)k+1 

Factorial series play a similar role in the theory of difference equations as inverse power series in the theory of 
differential equations, and classic books on finite difference such as the ones by Milne-Thomson |29], or Norlund 
032143411 ) treat factorial series. A contemporary discussion of the use of factorial series in the context of difference 



equations can be found in a recent article by Olde Daalhuis 113511 . 

But I am are much more interested in the convergence properties of factorial series, which fortunately can be 
analyzed easily. If we use [J, Eq. (6.L47) on p. 257] 

r(z + fl)/r(z + /7) = z"-''[l+0(l/z)], z^oo, (3.2) 

we obtain the asymptotic estimate n!/(z)„+i = 0(n^') as n — > oo. Thus, the factorial series ( ILlb converges with the 
possible exception of the points z — —m with m G No if and only if the associated Dirichlet series £l{z) — YZ=\ o-n/n^ 
converges (see for example [22, p. 262] or [23, p. 167]). Accordingly, a factorial series converges for sufficiently large 
5R(z) even if the reduced series coefficients a„ in (ILll i grow like a fixed power «" with a > as n ^^ oo. 

Factorial series are closely related to the beta function, which is usually defined as the following ratio of gamma 
functions [1, Eq. (6.2.2)]: 

^(^'^) = W^' ^'^e<^- ^3.3) 

Thus, the ratio n\/{z)n+\ can be expressed as a beta function: 

B{z,n + \) = j^. (3.4) 



Accordingly, a factorial series can also be expressed as an expansion in terms of beta functions 11291 p. 288]: 

n(z) = Y^a,B{z,n + \). (3.5) 

n=0 

The beta function possesses numerous integral representations. For our purposes the most useful one is the so- 
called Euler integral of the first kind (see for example [i Eq. (6.2. 1)]): 

B{x,y) = f t'^-^il-ty-^At, 5R(x),9^(y)>0. (3.6) 

4 



Combination of (13.5b and ( 13.61 ) yields the following integral representation: 

n! r' 



{z)n+\ Jo 



t'-'{l-tydt, 5R(z)>0 neNo. (3.7) 



If we now combine ( 13.4b . ( 13.5b and ( 13.7b and interchange integration and summation, we obtain the following integral 
representation (see for instance |31, Satz I on p. 244] or |29, p. 289]): 



n(z) = / r-^(p{t)dt, 3i{z)>Q, 
Jo 



(3.8a) 



(Pit) =. ^fl„(l-f)". (3.8b) 

This integral representation is of considerable importance: Frequently, the properties of Q.{z) can be studied more 
conveniently via the corresponding properties of the conjugated function (p{t) than via the defining factorial series 
(II. lb (see for example [31, Kapitel XVII]). As discussed in Section[8] this integral representation can also be used for 
the evaluation of factorial series. 

4. Transformations of inverse power series and factorial series 

Inverse powers l/z*^+' and inverse Pochhammer symbols l/{z)k+i can be transformed into each other via (lA.llb 
and ( IA.14b . respectively. Therefore, inverse power series and factorial series can also be transformed into each other. 
This has been known for a very long time. The algebraic processes effecting these transformations were already 
described in Nielsen's book Bill which was first published in 1906. It seems, however, that these potentially very 
useful transformation formulas are now largely forgotten. 

Let us assume that a function y/: C — s- C possesses the inverse power series ^/(z) = L"=o<^n/z"^'- If we insert 
(IA.14b into this series and rearrange the order of summations, we obtain the transformation formula 

Li^-Efer7L(-i)''s(^'(-'M)-M, (4.1) 

n=0 ^ m=0 K^/m+l ^=o 

which shows that the coefficients of the factorial series are weighted averages of the power series coefficients involving 



Stirling numbers of the first kind. I employed this transformation already in 115311 in order to speed up the convergence 
of truncated asymptotic expansions for the truncation errors of series expansions for special functions. 

An inverse expansion can also be derived. Let us assume that a function X'- 'C ^ C possesses a factorial series 
X{z) = lln=Q'^"/i^)n+i- If ^^ insert ( lA.llb into this series and rearrange the order of summations, we obtain the 
transformation formula 

L7^-L^-^Li-^r^''H'n,H)d,. (4.2) 

„=0 V^/n+l m=0 ^ 11=0 

The operations producing the transformation formulas ( 14.1b and (14.2b are purely formal. Therefore, we cannot 
tacitly assume that the inverse power or the factorial series necessarily converge. This has to be checked explicitly in 
each case. 

The sign patterns in the inner sums ^^o(-l)''S('^(m,^)c^ and ^'^o(-l)''S(^^(m,;^)c/^ in (gJj and ^^ are of 
crucial importance for the convergence or divergence of the formal expansions ( 14.1b and ( 14.2b . If the signs of the terms 
in the inner sums alternate, we can hope for a substantial cancellation as in binomial sums L/=o(~l)^( )/«+;' but if 
all terms have the same sign, a potentially explosive accumulation can take place. In the former case, convergence is 
likely, while in the latter case we should be prepared for divergence. 

The ability of the Stirling numbers of the first kind to achieve a cancellation is immediately obvious from its finite 
generating function. Setting z — k with 1 < A: < n — 1 in (lA.lb yields Lv=o^^S ("j v) — for n > 2. Let us now 
assume that the coefficients c„ in ( 14.1b have strictly alternating signs. Consequently, ( IA.3b implies that the terms of the 
inner sum (— 1)'"^'" o(— l)''S''-'(OT,/i)c^ — Lu=o 1^ (w,M)|c^ in (14.1b also have strictly alternating signs, and we 
can hope for a substantial cancellation. The inner sum should for large values of the outer index m be much smaller in 



magnitude than its individual terms. Because of this cancellation, we can hope that a function defined by a divergent 
inverse power series with strictly alternating coefficients can be expressed and computed by a convergent factorial 
series. 

The situation is not nearly as nice if the coefficients c„ in ( 14.11 ) all have the same sign. If we set in ( IA.2I ) z — 
I and use (lAJt . we obtain (-l)'"^'^o(-l)''S(')(m,^) = I:™=o|S'''('«,M)| = ml. Accordingly, the coefficients 
T!Il=o\^ (m, IJ-)\c n of the factorial series should be (much) larger than the coefficients c„ of the inverse power series. 

5. Stieltjes series 

A function F : C — J- C is called a Stieltjes function if it can be expressed by the Stieltjes integral 

^W = / -^^ |arg(z)|<;r. (5.1) 



Here, <P{t) is a bounded, nondecreasing function taking infinitely many different values on the interval < f < oo. 
Moreover, the moment integrals 



H„ = / t"d^{t), neNo, (5.2) 

Jo 

must be positive and finite for all finite values of n. 

Detailed discussions of Stieltjes series and their special role in the theory of summability can be found in the books 
by Bender and Orszag [3., Chapter 8.6] or Baker and Graves-Morris [2, Chapter 5]. In the case of divergent Stieltjes 
series, it can be shown rigorously that the Pade approximants [n + j/n] with fixed y > — 1 converge to a uniquely 
determined Stieltjes function as n — > oo. 

An inverse power series representation for such a Stieltjes function F{z) can be derived easily. We insert the 
geometric series Lv=o(~0^/-Z*'^' — ^/{z + t), which converges for \t/z\ < 1, into the integral representation (15.1b and 
- ignoring all questions of legitimacy and convergence - integrate term-wise from to oo using ( 15.2b . Thus, a Stieltjes 
function F(z) can at least formally be represented by its Stieltjes series 

^W-L^^- (5-3) 



v=0 



z" 



A factorial series for a a Stieltjes function F{z) can also be derived quite easily. For that purpose, we use the 
convergent factorial series (IA.12b with w = — f in the integral representation (15.1b and interchange integration and 
summation: 

nz)^ riT^dOW (5.4) 

^ ' {-t)nd^it). (5.5) 



If we now expand (— f)„ via ( lA.lb and do the resulting moment integrals according to (15.2b . we obtain: 

^W = E 7T^ E S<''(«,v) rt-d<P{t) (5.6) 

„=0 \Z)n+l v=0 •'0 

= E7^Es"'(«,v)Mv. (5.7) 

„=0 \Z)n+l v=0 



Obviously, the factorial series (15.71 ) for a Stieltjes function is a special case of the more general result (14.1b . 

Again, cancellation is the reason why we can hope that the factorial series (15.7b converges and is computationally 
useful even if the corresponding Stieltjes series ( 15.3b diverges. The power t" in the moment integral ( 15.2b is positive, 
but the Pochhammer symbol (— f)„ in the generalized moment integrals /J°(— f)„d<I>(f) in (15.51 ) has zeros for f = 
0, 1 , ... n — 1 and alternates in sign. This leads to a substantial cancellation. 



6. Summation of the asymptotic series for the exponential integral 

The practical usefulness of the transformation of a factorially divergent inverse power series to a convergent 
factorial series via (14.11 ) or via ( 15.71 ) in the case of a Stieltjes series can be demonstrated by means of the exponential 
integral II Eq. (5.1.1)] 

/■°° exp(-f)df 
Ei{z) = I -^V^' ^ ^^ 

n 

which possesses the following asymptotic expansion as z — )► oo [ 1 , Eq. (5.1.51)], 

ze-^£i(z) ~ ^ ^4^ = 2^o(l,l;-lA), z^-- (6.2) 

m=0 ^ 

The exponential integral can also be expressed as a Stieltjes integral yj, Eq. (5.1.28)]: 

eE.iz) = /" ^^^P^^^ = i r ^^^^H^ . (6.3) 

^ ^ Jo z + f z Jo l+t/z 

Thus, the divergent inverse power series ( I6.2l i is a Stieltjes series. 
If we now combine (16.21 ) with ( 15.61 ) and ( 15.71 ). we obtain: 



„=0 [Z)i'+l Jo „=o UJ«+1 v=0 



(6.4) 



As discussed in Sections|4]and|5] the transformation of an inverse power series to a factorial series only produces 
a numerically useful result if a substantial cancellation takes place in the inner sums in ( 14.11 ) or (15.7b . For that purpose, 
the inverse power series coefficients (— l)"n! and the factorial series coefficients (— l)"Lv=o^ ("'^)^' are displayed 
in Table [U 

Table 1 : Leading coefficients of the inverse power and the factorial series of the exponential integral. 

~o i r 

1 -1 -1 

2 2 1 

3 -6 -2 

4 24 4 

5 -120 -14 

6 720 38 

7 -5 040 -216 

8 40 320 600 

9 -362 880 -6 240 

10 3 628 800 9 552 

11 -39 916 800 -319 296 

12 479 001600 -519 312 

13 -6 227 020 800 -28 108 560 

14 87 178 291200 -176 474 352 



The numbers displayed in Table[T]show that there is indeed a substantial amount of cancellation in the inner sums 
in ( 16.4b : The coefficients of the inverse power series grow much more rapidly in magnitude than the coefficients of the 



factorial series. This is enough to produce a convergent resuh. For example, the first 15 coefficients of the factorial 
series produce for z = 5 the following result: 

14 /_! y; n 

Ey7Es(')(«,v)v! 

"=" ^ ' ;=° , , = 1.000 000 764 (6.5) 

exp(5)£'i(5) 

7. Conversion of a power series to a factorial series 

Assume that a function /: C — s> C possesses a power series f{z) = L,T=o 7nZ". For a transformation of this power 
series to a factorial series, we express it as an inverse power series in 1 /z: 



1£ 7. 
«=o 

If we now use (14. Il l, we obtain a factorial series in 1 /z: 



■^^^^ Z ?n (1/Z)«+1 



(7.1) 



L r«2" = - E 777^ E (-irs(')KM)7,. (7.2) 

„=0 Z ,„^o l^i/Zjm+l ^^0 

An equivalent factorial series was considered by Ramis and Thomann [i36i p. 20]. Thomann (|42, p. 526] and 1431 
Section 5.3]) considered similar expansions in terms of generalized factorial series, where the Pochhammer symbols 
are replaced by products z(z + fi ) ... (z + f„). It is clear that generalized factorial series are at least potentially more 
powerful than their ordinary counterparts. However, it is not a priori clear how the parameters {tn]n=i should be cho- 
sen. In addition, Thomann's formulas contain so-called generalized Stirling numbers instead of the ordinary Stirling 
numbers of the first kind. Therefore, it is not immediately obvious whether Thomann's generaUzed transformation 
formula is really more useful than ( 17.2b . 

Further manipulations of the Pochhammer symbol in (I7.2l i are possible: 

(7.3) 



1 


1 z^+1 z f. z 


(lA)m+l 


nu[k+ 1 A] nr=i k[z+i/k\ m\ |-J ,+ i/k 


Inserting this into \l .2\ yields: 




Ervz^ 

v=0 


oo / 1 yji m 111 

-E^n^E(-irs('W)7.. 

,„=0 "^- k=l ^+ ^/'^ Li=Q 



(7.4) 

If the power series coefficients jv have strictly alternating signs, we can expect cancellation in the inner sum involving 
the Stirling numbers. Moreover, for z > we have z/{z + l/k) < 1 for A: G N. Since the transformation ( 17.3b also 
produces the factorial I /ml, we can expect that the transformation ( 17.4b produces a convergent and numerically useful 
result even if the coefficients jv diverge factorially in magnitude. 

8. Tlie quartic anliarmonic oscillator 

In their seminal articles ||4J,|5|], Bender and Wu showed that the Rayleigh-Schrodinger perturbation expansions for 
the energy eigenvalues E^'">{^) of the anharmonic oscillators defined by the Hamiltonians 

H^'"\li) ^ p^+.f + pf"', m = 2,3,4,..., P^-i-^, (8.1) 

ax 

diverge quite violently for every nonzero coupling constant /3 (here, the same notation as in i50ll5Ul55ll is used). Later, 



perturbation expansions with a similar type of divergence were discovered in the case of other quantum mechanical 



systems, which ultimately created a new sub-discipUne of theoretical physics called large order perturbation theory 
(see for example the book by Le Guillou and Zinn-Justin 124] and the articles reprinted there). 

In the following years, a lot of work has been done on the summation of divergent perturbation expansions (a 
mathematically oriented overview can be found in OSll ). In particular the quartic anharmonic oscillator with m = 2 
in ( 18. lb has become a very popular computational laboratory for theoretical physicists. In spite of its simplicity, 
the quartic anharmonic oscillator leads to challenging computational and conceptual problems, as documented in 
countless articles (far too many to be cited here). 

In this article, I am exclusively interested in the summation of the divergent Rayleigh-Schrodinger perturbation 
expansion 

£(2)(^) = £/,(2)^« (8.2) 

«=0 

for the ground state energy of the quartic anharmonic oscillator by transforming it to a factorial series via ( I7.4l i. Long, 
but nevertheless incomplete lists of references dealing with other approaches for the summation of the divergent 
perturbation expansions of the anharmonic oscillators can be found in OSOIdU 15511 or also in 11491 Kap. 10]. 



If the convention ( 18. lb for the Hamiltonian is used, the coefficients b„ possess in the case of large indices n the 
following leading order asymptotics (see for example 115 OL Eq. (2.3)]): 

Z,(2) ^ („l)«+i(|^r(« + l/2)(3/2)«, «^-. (8.3) 

This asymptotic estimate shows that the perturbation series ( 18.2b diverges for all )3 7^ like the generalized hyperge- 
ometric series 2/^0(1/2, 1; -3i3/2) = i:^o(l/2)m(-3j3/2)"' |47, Eq. (1.10)]. 

It was shown rigorously bu Simon ll37l Theorem IV.2.1] that the perturbation expansion 

AE^^KP)^ib^'!iP" (8.4) 

11=0 

for the energy shift defined by E^^^P) = fe[)^' + j3A£'(2)(j3) = 1 + j3A£'P)(j3) is a Stieltjes series. This is a highly 
advantageous feature. As discussed in Section |5] this implies that the perturbation series (18.4b corresponds to a 
uniquely defined Stieltjes function since it is Pade summable. Moreover, the asymptotic estimat e (|8.3b implies that 
the terms of the perturbation expansion ( 18.41 ) have for )3 > strictly alternating signs (see also 1491 Tabelle 10-1]), 
which is advantageous if we want to sum it with the help of (17.4b . 

If we transform the perturbation series ( 18.41 ) for the energy shift with the help of ( 17.4b . we obtain the following 
expansion for the ground state energy of the quartic oscillator: 

E^^\P) = 1+/3 i ^-^Ujlmti-^rS^'H'n,^)b%. (8.5) 

m=0 "*• k=l f^^ ^/'^ 11=0 

(2) 

The first 34 terms of the infinite series on the right-hand side yield for j3 = 1 /5 the energy Ep^ (1/5) = 1.118305..., 

(2) 
which is less accurate than the energy £'p^( 1/5) = 1.118 292 654 373... obtained by computing the Pade approx- 

imants [17/17] from the first 34 terms of the perturbation expansion ( 18.4b for the energy shift A£''^'(l/5). These 

(2) 
approximations can be compared to the "exact" energy Zigxactl 1/5) = 1.118 292 654 367 039 154... obtained by a 

very sophisticated summation calculation [49, Tabelle 10-9]. 

Thus, the truncated series expansion ( 18.51 ). which does a linear transformation of the perturbation series coeffi- 
cients b\,b2, ■■■, ^34, is less efficient than the highly nonlinear Pade approximant using the same number of coeffi- 
cients b^ . Nevertheless, improvements are possible if we use the integral representation ( 13.8b for the evaluation of 
the corresponding factorial series i2(z) = Ljr=o'^«"V(-z)n+i- 

A direct use of a truncation of the power series (p{t) = Lr=o'^«(l ^ 0" defined by ( I3.8bb in the integral repre- 
sentation ( I3.8ab does not lead to an improvement since integration is linear It is, however, possible to replace the 
truncated power series for (p{t) by a Pade approximant in 1 — f and to evaluate the resulting expression by numerical 
quadrature. The use of the Pade approximant [17/17] to (p{t) in dlSl l yields £'intFs(l/5) = 1.118 292 654 369..., 



which is better than the direct Pade summation result £'p^(l/5) = 1.118 292 654 373..., but less accurate than 
£^Bp(l/5) = 11 18 292 654 367 039 152... obtained by doing a so-called Borel-Pade transformation originally 
introduced in |2l | . 

9. Outlook 

A single article cannot provide an exhaustive treatment of the numerical utilization of factorial series for the 
summation of divergent (inverse) power series. Because of space limitations, many interesting or potentially useful 
aspect of the theory of factorial series were treated only superficially or even completely ignored. Nevertheless, I hope 
that this article will inspire others. 

Inner sums of the type of (— 1 )"'£"()(— 1)''S' ■*(»?, /i)c^ occurring on the right-hand side of ( 14. Il l are the key 
quantities of this article. Recurrence formulas or even alternative closed form expressions would obviously be de- 
sirable. In the case of Stieltjes series, a recursive scheme for the computation of the generalized moment integrals 
./o° (^Ond4>(f) in ( 15.51 ) can be derived. 

The numerical examples of this article also raise questions. It is an obvious question whether and how well the 
convergence of factorial series can be accelerated by sequence transformations. When I looked at the convergence 
of the factorial series (16.41 1 for Ei (z), I applied sequence transformations to speed up its convergence. However, 
convergence was not improved substantially by the transformations I used. At the moment, it is unclear whether this 
is a specific feature of the factorial series ( I6.4l i. or whether we face a problem of a more general nature. It could 
be that the convergence of factorial series can only be accelerated effectively if other, specially designed sequence 
transformations are used. This should be investigated. 

It should be worthwhile to investigate whether the transformation formula (17.4b can also be used profitably in 
the case of convergent power series as a convergence acceleration tool. One can also hope that in the case of suffi- 
ciently simple power series coefficients 7„ explicit expressions for the inner sum (— 1)™52"'^q(— 1)''S' '(m,/i)7^ can 
be found, which would yield new explicit expressions in terms of factorial series for functions defined by power series. 

In Section[8] the truncated power series (p{t) = Lr=o'^n(l ~0" defined by ( I3.8bl ) was converted to a Fade approxi- 
mant in 1 — f , which was inserted into the integral representation (13.81 ). One should investigate whether other sequence 
transformations produce better approximations to (p{t) than Fade approximants. 

Appendix A. Stirling numbers 

With respect to notation, the theory of Stirling numbers is a mess. This is partly due to the fact that numerous 
different symbols are used in the literature (a discussion of the various notations can be found in fl, p. 822]). To 
make things worse, different and incompatible notations are used for factorial expressions. In special function theory, 
Fochhammer symbols (z)„ = z{z +l)...(z + n — l)=r(z + n)/r{z) are consistently used, but in combinatorics, it 
is more common to use instead falling factorials z(z — 1) . . . (z — « + 1) = r(z + l)/r(z — n + 1). Unfortunately, in 
the literature of combinatorics falling factorials are often denoted by the symbol (z)„ normally reserved in special 
function theory for Fochhammer symbols. 

The Stirling numbers S' ' (n, v) of the first kind are the polynomial coefficients of a Fochhammer symbol {z~n + 
1)„ = z(z - 1) . . . (z - n + 1) = r(z+ l)/r(z - n + 1) (see for example [4L Eq. (1) on p. 56]); 



(z-n+l)„ = (-l)"(-z)„ = ^S(')(n,v)^^ hgNq. (A.l) 

v=0 

If we use (z — « + 1 )„ = (— 1 )" (— z)„ and replace z by — z, we obtain: 

W« = (-1)" E (-l)^S(')(n,v)z^ «eNo. (A.2) 

v=0 

If z > holds, the coefficients of all powers z^ with < V < n in the expansion of {z)n are either zero or positive. 
Thus, ( IA.2I 1 implies 

(-l)«-^S(')(n,v) = |S(i)(«,v)|. (A.3) 
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The Stirling numbers S^ '(n, v) of the second kind are usually defined as follows (see for example 14 li Eq. (14) 
on p. 58]): 

z" = ^ S(2)(„,v)(z-V+l)v, neNo. (A.4) 

v=0 

The Stirling numbers of the first and second kind correspond to triangular matrices that transform the polynomial 
sequences {{z — n+l )«}"=o and {z" }~^o into each other Since these transformations must be invertible, orthogonality 
relationships exist. 

If we replace in ( lA.lb the powers on the right-hand side by Pochhammer symbols according to ( IA.4l l. we obtain: 



(z-«+l)„ - ^S(')(«,v)£sP)(v,^)(z-^+l),, „eNo. (A.5) 

V=0 A:=0 

By interchanging the order of the summations summations, this expression can be rewritten as follows: 

(z-« + l)„ = X;(z-A:+1), ^S(')(«,v)S(2)(v,^), neNo. (A.6) 

A:=0 v=k 

Thus, we obtain the following well-known orthogonality relationship: 

f S(')(«,v)SP)(v,^) = 5„,, k,neNo. (A.7) 

v=k 



By replacing in the Pochhammer symbols on the right-hand side of ( IA.4I I by powers according to ( lA.ll i. we obtain 
the following alternative orthogonality relationship: 

'£S^^\n,v)S^'\v,k) = 5„k, k,neNo. (A.8) 

V=k 

The Stirling numbers S' '(n, v) of the second kind possess the following infinite generating function (see for 
example |41, Eq. (16) on p. 58]): 

1 °° 

il-t)il-2r)...il-kr) -]^f'^' + '^'^^'^ .eN, \t\<l/k. (A.9) 



The substitution t = 1/z yields: 

1 



(z-k) 



k+l K=0 



Y,S^^Hk+K,k)z-''-''-\ \z\>k. (A.IO) 



If we now use {z — k)ic+i = {—^Y^^{—z)k+i and replace z by — z, we obtain the following inverse power series 
expansion of an inverse Pochhammer symbol [3li Eq. (9) on p. 68]: 



A convenient starting point for the derivation of a factorial series for an inverse power l/z^+' is the following 
factorial series Oil Eq. (3) on p. 77]: 

^=£7^, 3iiz-w)>0. (A.12) 

Next, we apply d^[l/(z - w)]/dw^' = kl/{z- w)''+^ with k e No to (IA.12b . which yields: 



(z-w)*+i dw\toW«+i 
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The Pochhammer symbol (w)„ can be expanded with the help of ( lA.lb . Then, the differentiations can be done in 

i I 
closed form and we obtain - after setting w = - the following factorial series for an inverse power [31, Eq. (6) on p. 

78]: 

Appendix B. Sequence inversion by triangular orthogonal matrices 

Let us assume that we have two matrices A ~ {A,„„}„ „>o and B ~ {Bmn}m,n>Q, which are triangular, 

Anm = B,„n = 0, n>m, (B.l) 

and which satisfy the orthogonality relationships 

Y^B„rA,, = 5„k, (B.2) 

r=k 
n 

2_^A„rBi-k — dnk- (B.3) 

The orthogonality relationships ( IA.7l i and ( IA.8I ) involving the Stirling numbers of the first and second kind are an 
example of such a pair of orthogonality relationships. 

We also assume that there are two sequences {xn}'^=Q and {}'n}~=o whose elements are connected by finite linear 
combinations 

r 

yr = Y^ArkXk, re No, (B.4) 

k=Q 

and two sequences {Mn}"=o ^^^ {*^''"}r=0' whose elements are connected by infinite series expansions 

"r = Y^AnrWn-, ^ £ Nq . (B.5) 



In the book of Charalambides lllSL Example 8.2 on pp. 284 - 285] it is shown that the finite linear combination 



( IB.4l i then possesses the inverse relation 



Xn ^ Y. ^"kyk I n e No , (B.6) 

*:=0 



and that the infinite series expansion ( IB.5I I then possesses the inverse relation 

Wk - f^B^kUr, ^GNo. (B.7) 

r=k 

Obviously, the two finite linear combinations ( IB.4b and ( IB.6I ) generalize the finite generating functions ( lA.lb and 
( IA.4b of the Stirling numbers, and the two infinite series expansions dB.Sb and (IB. 7b generalize the series expansions 
( lA.llb and ( IA.14b connecting inverse powers and inverse Pochhammer symbols. 

Further generalizations are possible. Let us now consider the following, in general formal infinite series expansion: 

F = £rj„3;„. (B.8) 

n=0 

With the help of ( IB.4b . the elements of the sequence {ynj^^o '-^'^ ^^ replaced by the elements of the sequence {-^hI^^o' 

F = Y.ri„f^A„kXk. (B.9) 

«=0 k=0 
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The order of the summations of this expansion can be rearranged: 

F = Xo[77oAoo + J?lAio + T72A20 + ...] + Xy [7?lAii+772A2i+773A3i + ...] 

+ X2[tJ2A22 + T73A32 + T74A42 + ---] + -1^3 [t?3A33 + r]4A43 + ^5^53 + •■■ ] + ••• (B.IO) 

^ Y^x„ Y^Ak„i]k- (B.ll) 

n=0 k=n 

Thus, the substitution {yn]'^=Q -^ {-''^n},7=o ^^ dB.Sb produces an ahernative expansion for F having the same general 
structure: 

F = Y, ^nXn, ^n = J^ Aa.„ TJ^. = ^ A„+v,„ TJ„+v . (B.12) 

j;=0 k=n v=0 



Unknowingly, I had used transformation formulas for expansion coefficients of the kind of ( IB.12I I already in 154 1 
when I studied the transformation of Laguerre series to power series. 

Next, we consider the following, in general formal infinite series expansion: 

G = £ V„u„. (B.13) 

«=o 

With the help of ( IB. 5b . the elements of {««}"=() can be replaced by the elements of {vv„}~^q: 

G = £u„£A,„Wi. (B.14) 

«=0 k=n 

The order of the summations of this expansion can be rearranged: 

G = woU„Aoo + vv'i [AioUo+AiiUi] + W2 [A20U0+A21 Ui +A22U2] +■■■ (B.15) 

= L ^" Jl^^kVk- (B.16) 

n=0 k=0 

Thus, the substitution {m„}~^o ^ {'^«}r=o ^^ < IB.13l l produces an alternative expansion of the same general structure: 

CO u 

G = Y, (^'^n, co„ = ^A„i,%. (B.17) 

11=0 k=0 

These relationships for G obviously generalize the transformation formulas (14.11) and (14.21 ). which transform factorial 
series and inverse power series into each other 

Let me emphasize once more that all operations considered in this Section are purely formal. Accordingly, we 
cannot tacitly assume that the transformation formulas necessarily lead to convergent expansions. This has to be 
checked explicitly in each case. 

This Appendix was inspired by the orthogonality relationships ( IA.7b and ( lA.Sb involving the Stirling numbers 
of the first and second kind. But many other mathematical objects possess similar features. Recently, I had studied 
the transformation of Laguerre expansions to power series expansions [54]. As is well known, generalized Laguerre 
polynomials and powers are connected by finite sums of the type of ( IB. 4b and (IB. 6b which can be inverted. Accord- 
ingly, these transformation formulas possess certain orthogonality properties when written in matrix form. Moreover, 
the coefficients of Laguerre series and of power series are connected by infinite series expansions of the type of 
(IB. 121 ). I only understood the wider significance of these features after 1 54] was published when I studied the book by 
Charalambides llisll more carefully. 

[1] M. Abramowitz, I. A. Stegun (eds.). Handbook of Mathematical Functions, National Bureau of Standards, Washington, D. C, 1972. 

[2] G. A. Baker, Jr., P. Graves-Moms, Pade Approximants, 2nd ed., Cambridge U. P., Cambridge, 1996. 

[3] C. M. Bender, S. A. Orszag, Advanced Mathematical Methods for Scientists and Engineers, McGraw-Hill, New York, 1978. 

[4] C. M. Bender, T. T. Wu, Anharmonic oscillator, Phys. Rev. 184 (1969) 1231 - 1260. 

[5] C. M. Bender, T. T. Wu, Large-order behavior of perturbation theory, Phys. Rev. Lett. 27 (1971) 461 - 465. 

13 



[6] E. Borel, Lefons sur les Series Divergentes, 2nd ed.. Editions Jacques Gabay, Paris, 1988. Originally published by Gauder-Villars, Paris, 

1928. English translation by C. L. Critchfield and A. Vakar, Lectures on Divergent Series, Translation LA-6140-TR, Los Alamos Scientific 

Laboratory, Los Alamos, 1975. 
[7] R. Borghi, Evaluation of diffraction catastrophes by using Weniger transformation. Opt. Lett. 32 (2007) 226 - 228. 
[8] R. Borghi, Joint use of the Weniger transformation and hyperasymptotics for accurate asymptotic evaluations of a class of saddle-point 

integrals, Phys. Rev. E 78 (2008) 026703-1 - 026703-1 1. 
[9] R. Borghi, On the numerical evaluation of cuspoid diffraction catastrophes, J. Opt. Soc. Amer A 25 (2008) 1682 - 1690. 
[10] R. Borghi, Summing Pauli asymptotic series to solve the wedge problem, J. Opt. Soc. Amer. A 25 (2008) 211-218. 
[11] R. Borghi, Joint use of the Weniger transformation and hyperasymptotics for accurate asymptotic evaluations of a class of saddle-point 

integrals. II. Higher-order transformations, Phys. Rev. E 80 (2009) 016704-1 - 016704-15. 
[12] R. Borghi, M. Alonso, Efficient evaluation of far-field asymptotic series, in: ICTON '07: 9th International Conference on Transparent Optical 

Networks, vol. 3, IEEE Xplore, Piscataway, NJ, 2007. 
[13] R. Borghi, M. Santarsiero, Summing Lax series for nonparaxial beam propagation. Opt. Lett. 28 (2003) 774 - 776. 
[14] C. Brezinski, A. C. Matos, A derivation of extrapolation algorithms based on error estimates, J. Comput. Appl. Math. 66 (1996) 5 - 26. 
[15] C. Brezinski, M. Redivo Zaglia, A general extrapolation procedure revisited. Adv. Comput. Math. 2 (1994) 461 -477. 
[16] C. Brezinski, M. Redivo Zaglia, On the kernel of sequence transformations, Appl. Numer Math. 16 (1994) 239 — 244. 
[17] E. Caliceti, M. Meyer-Hermann, P. Ribeca, A. Surzhykov, U. D. Jentschura, From useful algorithms for slowly convergent series to physical 

predictions based on divergent perturbative expansions, Phys. Rep. 446 (2007) 1 - 96. 
[18] C. A. Charalambides, Enumerative Combinatorics, Chapman & Hall, Boca Raton, 2002. 
[19] J. Ci'zek, J. Zamastil, L. Skala, New summation technique for rapidly divergent perturbation series. Hydrogen atom in magnetic field, J. Math. 

Phys. 44 (2003) 962 - 968. 
[20] E. Delabaere, J.-M. Rasoamanana, Sommation effective d'une somme Borel par series de factorielles, Annal. I'lnst. Fourier 57 (2007) 421 - 

456. 
[21] S. Graffi, V. Grecchi, B. Simon, Borel summability: Application to the anharmonic oscillator, Phys. Lett. B 32 (1970) 631 - 634. 
[22] K. Knopp, Theorie und Anwendung der unendlichen Reihen, Springer- Verlag, Berlin, 1964. 
[23] E. Landau, Uber die Grundlagen der Theorie der Fakultatenreihen, Sitzungsb. Konigl. Bay. Akad. Wissensch. Miinchen, math. -phys. Kl. 36 

(1906) 151-218. 
[24] J. C. Le Guillou, J. Zinn- Justin (eds.), Large-Order Behaviour of Perturbation Theory, North-Holland, Amsterdam, 1990. 
[25] D. Levin, Development of non-linear transformations for improving convergence of sequences. Int. J. Comput. Math. B 3 (1973) 371 - 388. 
[26] J. Li, W. Zang, J. Tian, Simulation of Gaussian laser beams and electron dynamics by Weniger transformation method. Opt. Expr. 17 (2009) 

4959-4969. 
[27] J.-X. Li, W. Zang, Y.-D. Li, J. Tian, Acceleration of electrons by a tightly focused intense laser beam. Opt. Expr 17 (2009) 11850 - 11859. 
[28] B. Malgrange, Sommation des sries divergentes. Expo. Math. 13 (1995) 163 - 222. 

[29] L. M. Milne-Thomson, The Calculus of Finite Differences, Chelsea, New York, 1981. Originally published by Macmillan, London, 1933. 
[30] F. Nevanlinna, Zur Theorie asymptotischer Potenzreihen, Ann. Acad. Sci. Fenn. Ser A 12 (1919) 1-81. 
[31] N. Nielsen, Die Gammafunktion, Chelsea, New York, 1965. Originally published by Teubner, Leipzig and Berlin, 1906. 
[32] N.-E. Norlund, Lemons sur les Series d'lnterpolation, Gautier-Villars, Paris, 1926. 
[33] N.-E. Norlund, Lemons sur les Equations Lineaires aux Differences Finies, Gautier-Villars, Paris, 1929. 

[34] N. E. Norlund, Vorlesungen liber Differenzenrechnung, Chelsea, New York, 1954. Originally published by Springer- Verlag, Berlin, 1924. 
[35] A. B. Olde Daalhuis, Inverse factorial-series solutions of difference equations, Proc. Edinb. Math. Soc. 47 (2004) 421 - 448. 
[36] J. P. Ramis, J. Thomann, Some comments about the numerical utilization of factorial series, in: J. Delia Dora, J. Demongeot, B. Lacolle 

(eds.), NumericaL Methods in the Study of Critical Phenomena, Springer- Verlag, Berlin, 1981, pp. 12 - 25. 
[37] B. Simon, Coupling constant analyticity for the anharmonic oscillator, Ann. Phys. (NY) 58 (1970) 76 - 136. 
[38] B. Simon, Fifty years of eigenvalue perturbation theory. Bull. Amer Math. Soc. 24 (1991) 303 - 319. 
[39] D. A. Smith, W. F. Ford, Acceleration of linear and logarithmic convergence, SIAM J. Numer Anal. 16 (1979) 223 - 240. 
[40] D. A. Smith, W. F. Ford, Numerical comparisons of nonlinear convergence accelerators. Math. Comput. 38 (1982) 481 - 499. 
[41] H. M. Srivastava, J. Choi, Series Associated with the Zeta and Related Functions, Kluwer, Dordrecht, 2001. 
[42] J. Thomann, Resommation des series formelles, Numer Math. 58 (1990) 503 - 535. 

[43] J. Thomann, Procedes formels et numeriques de sommation de series d'equations differentielles. Expo. Math. 13 (1995) 223 - 246. 
[44] I. Tweddle, James Stirling's Methodus Differentialis: An Annotated Translation of Stirling's Text, Springer- Verlag, London, 2003. 
[45] G. Watson, The transformation of an asymptotic series into a convergent series of inverse factorials. Rend. Circ. Mat. Palermo 34 (1912) 41 

-88. 
[46] E. J. Weniger, Nonlinear sequence transformations for the acceleration of convergence and the summation of divergent series, Comput. Phys. 

Rep. 10 (1989) 189-371, Los Alamos Preprint math-ph/0306302 (http://arXlv.org). 
[47] E. J. Weniger, On the summation of some divergent hypergeometric series and related perturbation expansions, J. Comput. Appl. Math. 32 

(1990) 291 - 300. 
[48] E. J. Weniger, Interpolation between sequence transformations, Numer Algor 3 (1992) 477 - 486. 
[49] E. J. Weniger, Verallgemeinerte Summationsprozesse als numerische Hilfsmittel fiir quantenmechanische und quantenchemische Rech- 

nungen, Habilitation thesis, Fachbereich Chemie und Pharmazie, Universitat Regensburg, Los Alamos Preprint math-ph/0306048 

(http://arXlv.org) (1994). 
[50] E. J. Weniger, A convergent renormalized strong coupling perturbation expansion for the ground state energy of the quartic, sextic, and octic 

anharmonic oscillator, Ann. Phys. (NY) 246 (1996) 133 - 165. 
[51] E. J. Weniger, Construction of the strong coupling expansion for the ground state energy of the quartic, sextic and octic anharmonic oscillator 

via a renormalized strong coupling expansion, Phys. Rev. Lett. 77 (1996) 2859 - 2862. 
[52] E. J. Weniger, Mathematical properties of a new Levin-type sequence transformation introduced by Cizek, Zamastil, and Skala. I. Algebraic 

14 



theory, J. Math. Phys. 45 (2004) 1209 - 1246. 
[53] E. J. Weniger, Asymptotic approximations to truncation errors of series representations for special functions, in: A. Iske, J. Levesley (eds.). 

Algorithms for Approximation, Springer- Verlag, Berlin, 2007. 
[54] E. J. Weniger, On the analyticity of Laguerre series, J. Phys. A 41 (2008) 425207-1 - 425207^3. 
[55] E. J. Weniger, J. Ci'zek, F. Vinette, The summation of the ordinary and renormalized perturbation series for the ground state energy of the 

quartic, sextic, and octic anharmonic oscillators using nonlinear sequence transformations, J. Math. Phys. 34 (1993) 571 - 609. 
[56] E. J. Weniger, E. O. Steinborn, Nonlinear sequence transformations for the efficient evaluation of auxiliary functions for GTO molecular 

integrals, in: M. Defranceschi, J. Delhalle (eds.). Numerical Determination of the Electronic Structure of Atoms, Diatomic and Polyatomic 

Molecules, Kluwer, Dordrecht, 1989, pp. 341 - 346. 



15 



